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We study some discrete symmetries of unbiased (Hadamard) and biased quantum walk on a 
line, which are shown to hold even when the quantum walker is subjected to environmental effects. 
The noise models considered in order to account for these effects are the phase flip, bit flip and 
generalized amplitude damping channels. The numerical solutions are obtained by evolving the 
density matrix, but the persistence of the symmetries in the presence of noise is proved using 
the quantum trajectories approach. We also briefly extend these studies to quantum walk on a 
cycle. These investigations can be relevant to the implementation of quantum walks in various 

■ known physical systems. We discuss the implementation in the case of NMR quantum information 
' processor and ultra cold atoms. 

o : 
3 : 

1—5 . I. INTRODUCTION 

o : 

■ Random walk, which has found applications in many fields [l|, Q , is an important constituent of information theory. 
It has played a very prominent role in classical computation. Markov chain simulation, which has emerged as a 
powerful algorithmic tool Q, as well as many other classical algorithms, are based on random walks. Quantum 
random walks 0], the generalization of classical random walks to situations where the quantum uncertainties play a 

qq , predominant role, are of both mathematical and experimental interest and have been investigated by a number of 
—i ■ groups. It is believed that exploring quantum random walks [5| allows, in a similar way, a search for new quantum 
f**«« \ algorithms, a few of which have already been proposed 0, 0,111, HI- Experimental implementation of the quantum 

■ walk has also been reporte d flfj. 
quantum walks [H El Q \mM I ■ 

The evolution of a discrete classical random walk, involving steps of a given length, is described in terms of 
probabilities. On the other hand, the evolution of a discrete quantum random walk is described in terms of probability 
Oh, amplitudes. An unbiased one dimensional classical random walk with the particle initially at xq evolves in such a 

■ way that at each step, the particle moves with probability 1/2 one step to the left or right. In a quantum mechanical 
analog the state of the particle evolves at each step into a coherent superposition of moving one step to the right and 
one step to the left. 

q-( In the one dimensional quantum (Hadamard) walk the particle initially prepared in a product state of the coin 
(internal) and position (external) degree of freedom, is subjected to a rotation in the coin Hilbert space, followed 
by a conditional shift operation, to evolve it into a superposition in the position space. The process is iterated 
without resorting to intermediate measurements to realize a large number of steps of quantum walk, before a final 
' measurement. 

A one dimensional quantum (Hadamard) walk starting from initial state |0) ® \xo) or |1) <E> \xq) (where the first 
register refers to the coin degree of freedom, and the second register to the external, spatial degree of freedom) induces 
an asymmetric probability distribution of finding the particle after N number of steps of quantum walk. A particle 
with initial coin state |0) drifts to the right (solid line in Figure [l} and particle with initial state |1) drifts to the left 
(dashed line in Figure [1]). This asymmetry arises from the fact that the Hadamard coin treats the two directions |0) 
and |1) differently; it multiplies the phase by -1 only in case of |0). It follows that to obtain left-right symmetry, one 
needs to start the coin in the state |0) + i\l) (Figure [2]). 

In this paper we report that the quantum walk- both the unbiased (Hadamard) as well as the biased- remains 
invariant under certain operations, i.e., the probability distribution of the walker's position remains unchanged under 
the inclusion of these operations at each step of the walk. We refer to these discrete operations as symmetries of the 
quantum walk. We further study these symmetries in the presence of environmental effects, modeled by various noise 
channels. These results were obtained by numerical integration rather than simulation. We extend these studies to 



ll| . Various other schemes have been proposed for the physical realization of the 
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FIG. 1: (color online) Probability distribution of the quantum walk with the initial state |0) (solid line) and initial state |1) 
(dashed line) on the position. The distribution is for 100 steps. 
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FIG. 2: (color online) The probability distribution of an unbiased walker, and a biased quantum walker using an SU(2) operator 
of the form ([6]) as the quantum coin toss, (a) Unbiased quantum walker (9 — 45°); (b) Biased quantum walker with 6 = 15°; 
(c) Biased quantum walker with 8 = 75°. The distribution is for 100 steps. 



quantum walk on a cycle, which can be conveniently generalized to more general graphs. It is shown that the above 
symmetries do not hold, in general, for a quantum walk on a cycle and hence for the closed graph, but leads to other 
interesting behavior. These observations can have important implications for a better insight into, and for simplifying 
certain implementations of, quantum walks. 

This paper is organized as follows. Section |TT] briefly recapitulates the theory of quantum random walk. Section ITlTl 
discusses quantum walk augmented by various symmetry operations, both in the case of a biased and unbiased quantum 
coin. The above observations are generalized to the case of a noisy quantum walk in Section IIV1 with Section IIV Al 
treating noise as phase-flip and bit-flip channels, and Section HV Bl treating noise as a generalized amplitude damping 
channel. Whereas the numerical results presented here, involve evolving the density operator of the system, we have 
found it convenient to explain the symmetries using quantum trajectories. In Section|V]we extend these studies to the 
quantum walk on the cycle which can be generalzed to all closed graphs in general. In Section I VII we show that the 
application of these ideas can help simplify certain experimental implementations of quantum walk. We also discuss 
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the physical systems, NMR quantum-information processor and ultracold atoms, where the results presented in this 
article can be applied. In Section [VIII we make our conclusions. 



II. QUANTUM RANDOM WALK 

Unlike classical random walk, two degrees of freedom are required for quantum random walks, the internal, coin 
degree of freedom, represented by the two-dimensional Hilbert space TL Cl spanned by the basis states |0) and |1), and 
the particle degree of freedom, represented by the position Hilbert space TL P , spanned by basis states \x), x € Z. The 
state of the total system is in the space TC = TL C <£> TL p . The internal state of the particle determines the direction of 
the particle's movement when the conditional unitary shift operator U is applied on the particle, whose initial state 
is given by a product state, for example, \^i n ) — 7§[|0) + ® |*x )> where 

U = |0><0|®J>-l><:r| + |l><l|® J>+l)<s| 

)i+|l)(l|®it. (!) 



Here A and A^ are unitary operators that are notationally reminiscent of annihilation and creation operations, 
respectively. The conditional shift can also be written as 

U = exp(-2iS z ®Pl), (2) 

P, being the momentum operator and S z , the operator corresponding to the step of length /. 

Conditioned on the internal state being |0) (|1)), the particle moves to the left (right), i.e., U\Q)® \x) = |0) ® \x — 1) 
and U\l) (g> \x) = |1) ® |a; + l). Application of U on spatially entangles the TL C and TL P and implements quantum 

(Hadamard) walk, 

U\* in ) = i=[|0) ® e- iPl + |1) ® e^ xa ). (3) 
Each step of the quantum (Hadamard) walk is composed of a Hadamard operation (rotation) H, 

H = ^(] \), (4) 



on the particle, bringing them to a superposition state with equal probability, such that, 

(H®l)\0,x) = -^[|0,a:) + |l,a:)], 

(H®l)\l,x) = -^[|0,a;)-|l,a:)], (5) 

and, a subsequent unitary conditional shift operation, U, which moves the particle into an entangled state in the 
position space. 

The probability amplitude distribution arising from the iterated application of W = U{H ® 1) is significantly 
different from the distribution of the classical walk after the first two steps 0) ■ If the coin initially is in a suitable 
superposition of |0) and |1) then the probability amplitude distribution after n steps of quantum walk will have 
two maxima symmetrically displaced from the starting point. The variance of quantum version is known to grow 
quadratically with number of steps n, a 2 oc n compared to the linear growth, a 2 oc n for the classical walk. 

Figure [2] shows the probability distribution of finding the position of the quantum walker using a biased walker 
with 9 = 15° and 8 = 75°, respectively, using a general SU(2) operator of the form, 

, . _ / e*« cos(0) e< sin(0) \ , , 

o, g - I g _ lC sin(0) cog(0) I , (b) 



where we have set £ = £ = 0. Note that H = £?(0,45°,0). Replacing the Hadamard coin with an arbitrary SU(2) 
operator yields a biased coin toss, with 8 < 45° or 8 > 45°. The value of a can be increased or decreased, respectively. 
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III. BIT FLIP AND PHASE FLIP SYMMETRIES IN QUANTUM WALK 

Consider the application of the modified conditional shift operator U' = |1)(0| ® A + |0) (1| <S> A* instead of Eq. (fTJ). 
In place of Eq. ([5]), one has: 

U' = (X®I)exp(-2iS z ®Pl), (7) 

where X is the Pauli X operator. Since U' = XU, i.e., it is equivalent to an application of bit flip following U, 
conditioned on the internal state being |0) (|1)), the particle will move to the left (right) and changes its internal state 
to |1) (|0)). Thus, U'\0) <g> |ar) = |1) ® \x - 1) and U'\l) <g> \x) = |0) <g> \x + 1). 

A relevant observation in this context is that there are physical systems where the implementation of U' is easier 
than that of U (la |. In that case, applying a compensatory bit flip on the internal state, after each application of 
U', reduces the modified quantum walk to the usual scheme. In all this would require (n — 1) compensatory bit flips, 
which adds to the complexity of the experimental realization. However, this additional complexity can be eliminated. 
For an unbiased quantum walk, applying a bit flip in each step can be shown to be equivalent to a spatial inversion of 
the position probability distribution. A quick way to see why bit flips (for an unbiased quantum walk) are harmless 
is to note that they are also equivalent to relabeling the edges of the graph on which the quantum walk takes place, 
so that each end of each edge has the same label [l7|. (It is worth noting that in Ref. [HI], bit flips are employed to 
improve the practical implementation of a quantum walk of atoms in an optical lattice.) We may in this sense call a 
bit flip together with spatial inversion a symmetry of the unbiased quantum walk on a line. To be specific, a quantum 
walk symmetry is any unitary operation which may uniformly augment each step of a quantum walk without affecting 
the position probability distribution. Experimentally, the symmetries are useful in identifying variants of a given 
quantum walk protocol that are equivalent to it. This motivates us to look for other (discrete) symmetries of the 
quantum walk, which we study below. We begin with Theorem [T] where we note four discrete symmetries, associated 
with the matrices B^ (j = 1,2,3,4) (cf. Eq. ©), of the quantum walk. Thereafter two of these symmetries, B^ 
and B^ 2 \ are identified with operations that are relevant from the perspective of physical implementation. It is an 
interesting open question of relevance to practical implementation of quantum walks, whether other such symmetries 
of the quantum walk exist. 

Theorem 1 If B in Eq. is replaced by any of B^\ B^ 2 \ B^ 3 \ or B^\ given by, 

(1) = / e l « cos(0) e* c sin(6>) \ (2) = / cos(0) sin(0) 

~ ^ e^-O sin(0) - e *<*-« cos(0) J ' ~\ e -^sm(6) -e l ^-«) cos(0) 

(3 ) = / e^+«) cos(0) e^+O s in(0) \ (4) = ( e«*+« cos(0) e< sin(0) \ - . 

~V e-^sin(0) -e- l «cos(0) ) ' ~ \ e 1 ^-^ sin{9) -e^ cos(0) J ' { ) 

the resulting position probability distribution of the quantum walk remains unchanged. 

Proof. With the notation B = {b jtk } and B&> = {6^}, we find bfj = b jk e ij<l> , b$ = b jk e lk(t >, bfj = b jk e f ^, 
b^J — bj k e 1 ^, where the matrix indices j,k take values and 1, i = and the overbar denotes a NOT 

operation (0 ^> 1). The state vector obtained, after n steps, using B and B^ as the coin rotation operations, are, 
respectively, given by 

|*0 = (UB) n \a,f3)= J2 b j«,jn- 1 ---b j2 , jl b jua \j n ^ + 2J-n), 
h,h,— tin 

|* 2 ) = (UB^r\a,0)= E ^-i"-6AJ 1 6ii,«(^) i "- 1+ "- +il+a |jn,i8 + 2J-n) J (9) 

h,h,— ,3n 

where J ~ ji + ■ ■ ■ + j n - Consider an arbitrary state \a, b) in the computational-and-position basis. Now, (a, b\^>i) = 
e lTl ^{a, 6|*2); where r\ = j n -i + ■ ■ ■ + .7i + which is fixed for given a and b, and determined by b = /3 + 2 J — n and 
j n = a. As a result, |(a, 6|*i) | 2 + |(a, 6|*i) | 2 = |(a, 6|* 2 ) | 2 + I (a, ^1*2) | 2 - A similar proof of invariance of the position 
distribution can be demonstrated to hold when B is replaced by one of the other B^'s (j = 2,3,4). On account of 
the linearity of quantum mechanics, the invariance of the walk statistics under exchange of the B^'s and B holds 
even when the initial state \ot,0) is replaced by a general superposition or a mixed state. I 
Interchanging B and the B^'s may be considered as a discrete symmetry operation G : B — > B* (where B* denotes 
any of the B^''s in Eq. ©), that leaves the positional probability distribution invariant. We express this by the 
statement that 

W~GW, (10) 
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where G refers to the application of G at each step of the walk, and W refers to the walk operation of evolving 
the initial state through n steps and then measuring in the position basis. Knowledge of this symmetry can help 
simplify practical quantum walks. Below we identify two of these quantum walk symmetry operations B «-> £?W and 
B <-> B^ 2 \ associated with physical operations of interest. 

We first consider the phase shift gate $(</>) = |0)(0| + e l ^|l)(l| as a symmetry operation of a quantum walk. In our 
model, the quantum operation for each step is augmented by the insertion of <&(</>) just after the operation UB. At 
each step, the walker evolves according to: 



U$ = 



1 

e^ 

1 





A 



A 




1 






it 



it 



cos(0) 
sin((9) 



i(4>-i) 



cos(0) 
sin(0) 



e iC - sin(0) 
-e~ l( - cos(0) 

e lC - sin(0) 
-e*W~« cos(0) 



(11) 



This is equivalent to replacing B by B^ x \ which, according to Theorem [TJ leaves the walker distribution invariant. 
Thus the operation $ (</>), applied at each step, is a symmetry of the quantum walk. 

As a special case, the phase flip operation (Z) applied at each step, obtained by setting (j> = tt, is a symmetry of 
the quantum walk. Representing the inclusion of operations $ or Z at each step of the walk by or Z, respectively, 
we express this symmetry by the statements: 



W ~ $W, 
W ~ ZW. 



(12a) 
(12b) 



Unlike the phase flip operation, bit flip is not a symmetry of the quantum walk on a line. However, the combined 
application of bit flip along with angular reflection R (0 — > ir/2 — 0, i.e., sin0 «-> cos0, and £ <-> — £) and parity P 
(i «-> it) turns out to be a symmetry operation. These three operations commute with each other. By the inclusion 
of PRX, the walker evolves by (\PRX]UB) n . At each step, the walker evolves according to: 



PR 



1 

1 






1 



i- 



A 



it 



it 

cos(0) 
sin(0) 



cos(0) 
sin(0) 

sin(0) 
cos(0) 



sin(0) 
cos(0) 



-e< 



(13) 



This is equivalent to replacing B by B^ with (j> = tt, which, according to Theorem [TJ should leave the walker 
distribution invariant. Thus the operation PRX applied at each step, is a symmetry of the quantum walk. It will be 
convenient henceforth to choose £ = ( = 0, so that R will simply correspond to the replacement 8 — > 7r/2 — 8. 

Representing the inclusion of operations P, R or X at each step of the walk by P, R or X, respectively, we express 
this symmetry by the statements: 



W 

xw 



PRXIU, 
PBW. 



(14a) 
(14b) 



Eq. (|14ap was proved immediately above. Eq. (|14b[) follows from Eq. (|14a|) . since the operations P, R and X mutually 
commute, and X 2 = I. It expresses that fact that applying the X operation at each step is equivalent to replacing 
a quantum walk by its angle-reflected, spatially inverted counterpart. The observation made at the beginning of this 
Section pertains to the special case of 8 = 45°. By a similar technique the following symmetries may be proved: 



X.W ~ XZW 
The first equivalence easily follows from Eq. (|12p. 



ZXW. 



(15) 



IV. ENVIRONMENTAL EFFECTS 



A quantum walk implemented on a quantum computer is inevitably affected by errors caused by noise due to 
the environment. We consider three physically relevant models of noise: a phase flip channel (which is equivalent 
to a phase damping or purely dephasing channel), a bit flip channel and a generalized amplitude damping channel 
(T > 0). In all cases, our numerical implementation of these channels evolves the density matrix employing the 
Kraus operator representation for them. However to explain symmetry effects, it is convenient to use the quantum 
trajectories approach, discussed below. 
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FIG. 3: (color online) The effect of environmental decoherence on the position probability distribution of a biased quantum 
walker subjected to a noisy channel. Coin bias is of the form Eq. with 6 = 60° . The noise is modelled as a phase flip (solid 
line) or bit flip (dashed line) channel, characterized by Eqs. (|18|) and (|23|) . respectively, at various noise levels p: (a) p = 0.005 
(b) p = 0.05; (c) p = 0.1; (d) p = 0.5, which corresponds to a fully classical random walk. Comparing Figure (d) with Figure 
IHd), we note that the distribution in the case of maximal bit flip noise is the same. 



A. Decoherence via phase damping and bit flip channels 



In studying the status of the walk symmetries in the presence of noise, it is advantageous to employ the quantum 
trajectories approach [l8j . This simplifies the description of an open quantum system in terms of a stochastically 
evolving pure state, which allows us to adapt the symmetry results for the pure states, given in the preceding section, 
to mixed states. 

We call the sequence of walk step operations X = (U B n )(U B n -i) ■ ■ ■ (UB\) a 'quantum trajectory'. (More precisely, 
a trajectory refers to the sequence of states produced by these operations, for which the above serves as a convenient 
representation.) If all the -Bj's are the same, then X is the usual 'homogeneous' quantum walk W. In general, 
the Bj's may be different SU(2) operators, corresponding to a varying bias in the coin degree of freedom. More 
generally, each step of the walk may include generalized measurements whose outcomes are known (Section IIVB I) . If 
each walk step in X is subjected to a fixed symmetry operation G, the result is a new quantum trajectory GX = 
(U B^)(U B^i) ■ ■ ■ (UB*). We have the following generalization of Theorem [T] to inhomogeneous quantum walks on 
a line. 
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FIG. 4: (color online) The effect of environmental decoherence on the position probability distribution of a biased quantum 
walker subjected to a noisy channel. Coin bias is of the form Eq. ((SJl with 6 = 30° . The noise is modelled as a phase flip (solid 
line) or bit flip (dashed line) channel, characterized by Eqs. (|18|) and (|23|) . respectively, at various noise levels p: (a) p — 0.005 
(b) p = 0.05; (c) p — 0.1; (d) p — 0.5, which corresponds to a fully classical random walk. 



Theorem 2 Given any quantum walk trajectory X — (UB n ) ■ ■ ■ (UBi), the symmetry Q holds, i.e., X ~ GX . If the 
operation $ is restricted to Z, then the symmetries hold even when some of the U 's are replaced by U> 's. 

Proof. In the proof of Theorem [T] we note that if, in each step of the walk, we alter the rotation B by 
the transformation G, the proof still goes through. That is, l^i) = (UB n )(UB n -i) ■ ■ ■ (UB{)\a,(3) and \^>2) = 
(UB*)(UB*_ 1 ) ■ ■ ■ (UB*)\a, (3) produce the same position distribution. 

Suppose that in some of the walk steps, U is replaced with W . In place of Eq. we have 

|*i) = {UB n )---{U^B j )---{UB 1 )\a,P) 

= b in,jn-i ■ ■ ■ b 32ji b ji,a\jn,f3 + 2(Ji - J 2 ) - (ni - n 2 )), (16a) 

31,32,- ■■ 

|* 2 ) = (UBM)---(UiBf)...(UB{V)\a,P) 

= E K^-i---K,3^A-^y n - 1+ '-- +n+a \3niP + %h- J 2 )-(rai-n 2 )), (16b) 

31,32,- ■■ ,jn 

where J\ = ^2 k jk for the n\ steps fc where operator U is used, and J 2 = ^2iji for the n 2 steps / where operator U' 
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FIG. 5: (color online) Variation of standard deviation with noise level, for both phase noise and bit flip noise. Solid line is for 
an unbiased walker. For biased walker (a) 9 = 30° and (b) 9 — 60°. In the classical limit of p = 0.5, the standard deviation 
converges to a fixed value for bit flip noise, irrespective of 9, but different for phase flip noise. The convergence happens because, 
at maximum bit flip noise (p = 0.5), the measurement outcome in the computational basis is completely randomized, and the 
presence or absence of bias is irrelevant. On the other hand, the non-convergence in the case of phase flip noise is due to the 
fact that the asymptotic mixed state obtained via a phase damping channel depends on the initial state parameter 9. 

is used. Here J — J\ + J2 and n = n± + Ti2- Observe that the exponent of (—1) is effectively evaluated in modulo-2 
arithmetic. We can thus replace J by Jj — J2 in the exponent. Following the argument in Theorem [TJ we find that 
(a,b\^i) = e ie (a,b\^ 2 ), where = Ji - J 2 - a + a. ■ 

As a corollory, the symmetries Z and PRX hold good because they reduce to special cases of G. A question of 
practical interest is whether PRX and Z arc symmetries of a noisy quantum walk. Suppose we are given a noise 
process Af in the Kraus representation 



P 



E E iP E h 
3=0 



E4^ 



1. 



(17) 



With the inclusion of noise, each step of the quantum walk becomes augmented to (HUBk), where II is a random 
variable that takes Kraus operator values Ej. Thus, j\f corresponds to a mixture of upto n m trajectories or 'unravelings' 
Xi = (H(l n )UB n ) ■ ■ ■ (H(li)UBi), each occuring with some probability pi, where J^iPi = 1- ^ ^ an d PRX arc 
symmetries of an unraveling AT/, then the operations A/ and HA; = (H(l n )HUB n ) ■ ■ ■ (H(li)HUBi), where H denotes 
Z or PRX, must yield the same position probability distribution. In the case of bit-flip and phase-flip channels, there 
is a representation in which the £j-'s are proportional to unitary operators. Further: 

Theorem 3 If trajectories A; are individually symmetric under operation G, then so is any noisy quantum walk 
represented by a collection {Xi,pi}. 

Proof. The state of the system obtained via Af is a linear combination (the average) of states obtained via the Xj's. 
Thus, the invariance of the A^'s under G implies the invariance of the former. ■ 

This result, together with those from the preceding Section, can now be easily shown to imply that the symmetry 
H is preserved in the case of phase-flip and bit-flip channels. 

Decoherence via a purely dephasing channel, without any loss of energy, can be modeled as a phase flip channel 

S(p) = (l-p)p+pZpZ. (18) 

An example of a physical process that realizes Eq. (fT8|) is a two-level system interacting with its bath via a quantum 
non-demolition (QND) interaction given by the Hamiltonian 



H = H S + J2 hu k b{b k + H S J2 9k{b k + b\) +H 2 S J2 

k 



A 

tkOk 



(19) 
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Here H$ is the system Hamiltonian and the second term on the RHS of the above equation is the free Hamiltonian 
of the environment, while the third term is the system-reservoir interaction Hamiltonian. The last term on the RHS 
of Eq. (fl9|) is a renormalization inducing 'counter term'. Since [Hs, Hsr] — 0, Eq. (fl9|) is of QND type. 

Following Ref. [2(| (apart from a change in notation which switches |0) < — ► |1)), taking into account the effect 
of the environment modeled as a thermal bath, the reduced dynamics of the system can be obtained, which can be 
described using Bloch vectors as follows. Its action on an initial state: 



„ - ,s(l+W0)» (MO)) | (on) 

p0 =\ ( ff+ (o)) l(i -Wo))) 1, {) 



is given in the interaction picture by: 



£{P0) - {(a + (0))e-^) i(l -MO))) J' (21) 

The initial state (|20p may be mixed. (The derivation of the superopertor £ in terms of environmental parameters for 
the pure state case, given explicitly in Ref. [2pj | , is directly generalized to the case of an arbitrary mixture of pure 
states, since the environmental parameters are assumed to be independent of the system's state.) 

Comparing Eq. (|21[) with Eq. (|18p allows us to relate the noise level p in terms of physical parameters. In particular: 

p=i(l-exp[-(M 2 7W])- (22) 

When 7(f) « (either because the coupling with the environment is very weak or the interaction time is short or 
the temperature is low), p«0, tending towards the noiseless case. On the other hand, under strong coupling, j(t) 
is arbitrarily large, and p — > 1/2, the maximally noisy limit. The result of implementing channel (fl8|l is to drive the 
position probability distribution towards a classical Gaussian pattern [ItJ ■ The effect of increasing phase noise in the 
presence of biased walk is depicted in Figure U2 for the case of 9 — 60°, and in Figure 2J for the case of 9 = 30°. The 
onset of classicality is observed in the Gaussianization of the probability distribution. This is reflected also in the fall 
of standard deviation, as shown in Figure [51 

Decoherence can also be introduced by another noise model, the bit flip channel [l9l |: 

£(p) = (l-p)p + pX P X. (23) 

As with the phase damping channel, the bit flip channel also drives the probability distribution towards a classical, 
Gaussian pattern, with increasing noise [l7l |. The effect of increasing bit flip noise in the presence of biased walk is 
depicted in Figures [3] and 0J Here again, the onset of classicality is observed in the Gaussianization of the probability 
distribution, as well as in the fall of standard deviation, as shown in Figure [5] 

A difference in the classical limit of these two noise processes, as observed in Figure EJ is that whereas the standard 
deviation (in fact, the distribution) is unique in the case of the bit flip channel irrespective of bias, in the case of 
phase flip noise, the classical limit distribution is bias dependent. This is because phase flip noise leads, in the Bloch 
sphere picture, to a coplanar evolution of states towards the a z axis. Thus all initial pure states corresponding to a 
fixed 9 evolve asymptotically to the same mixed state [TiJ [2(| ■ This also explains the contrasting behavior of bit flip 
and phase flip noise with respect to bias, as seen by comparing Figures ([3]) and ([¥]). 

Representing the walker distribution by its standard deviation a, we may describe symmetry by the ratio of a 
without the symmetry operation to a with the symmetry operation. Figure [5] depicts the symmetry operation X for 
various bit flip noise levels. The convergence of the curves representing various #'s is a consequence of the complete 
randomization of the measured bit outcome in the computational basis. This implies that although X is not a 
symmetry of biased walk, it does become one in the fully classical limit. On the other hand, the symmetries PRX 
remain unaffected by noise. We note that, since the quantum walk here is evolved from the symmetric state |0) +z|l), 
and the bit flip and phase flip noise are not partial to the state |0) or |1), this is equivalent to setting P to 1, which 
explains the fact that the distributions in Figures [3] and [?] are spatially symmetric. Thus, RX by itself becomes a 
symmetry operation, which is manifested in the fact that in Figure [6] the values of the curve for complementary angles 
are the inverse of each other. Figure [7J shows that for either of the two noises, Z is a walk symmetry. 

Figure [5J depicts the symmetry of the RX operation at all phase flip noise levels, as evident from that fact that the 
values of the curve for complementary angles are the inverse of each other. From Figures [H [7J and [3J we note that for 
the Hadamard walk, all three symmetry operations Z, X and R are individually preserved. This is expected because 
here P = 1 as stated earlier, R = 1 by definition, so that the symmetry of PRX implies X = 1. 

With the inclusion of bias and an initial arbitrary state, the full symmetries Z and PRX would be required, as 
proved by the following theorem. 
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FIG. 6: (color online) Variation of the ratio of standard deviation without any symmetry operation to the bit flip symmetry 
operation with increasing bit flip noise level. 
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FIG. 7: (color online) Variation of the ratio of standard deviation without any symmetry operation to the phase flip symmetry 
operation with increasing bit flip (phase flip) noise level. 
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FIG. 8: (color online) Variation of the ratio of standard deviation without any symmetry operation to the bit symmetry 
operation with increasing phase flip noise level. 
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FIG. 9: (color online) Amplitude damping channel acting on a Hadamard walker at temperature T = 0. The distribution 
corresponding to intermediate values of p clearly show the breakdown of the RX symmetry. However, the extended symmetry, 
PRX (where P stands for parity operation (spatial inversion)) holds good. This is observed at all temperatures, (a) Probability 
distribution of finding the particle undergoing unbiased quantum walk on which amplitude damping channel is acting. This 
shows that even at T = 0, for sufficiently high coupling, the distribution turns classical, (b) Amplitude damping with bit flip 
symmetry. 




FIG. 10: (color online) PRX symmetry seen to hold in biased quantum walk subjected to amplitude damping (T — 0). The 
two cases are spatial inversions of each other. This holds for a generalized amplitude damping at any temperature, (a) walker 
with = 30° and bit flip; (b) 9 = 60°. 



Theorem 4 The operations PRX and Z are symmetries for the phase-flip and bit-flip channels. 

Proof. We may look upon the phase flip channel (| 18[) as a probabilistic mixture (in the discretized walk model) of 
2™ quantum trajectories with II 6 {I, Z}. By virtue of Theorem [31 it suffices to show that any given unravelling is 
invariant under Z and PRX. Consider an unravelling X 1 = ■ ■ ■ (IUB) (ZUB) (IUB) ••• = ••• (UB)(ZUB)(UB) 
This is the same as: ■ ■ ■ (UB)(UB')(UB) ■ ■ ■ , where B' = ZB, noting that Z commutes with U. Now, 
TLX 1 = ■■■ (IZUB)(ZZUB)(IZUB) ■■■ = ... (ZUB)(ZUB')(ZUB) which, by Theorem^ is equivalent to X 1 . 
Now, PRXI 1 = ■■■(PRX UB)(Z PRX UB)(PRX UB)--- = ■■■(PRX UB)(PR ZX UB)(PRX UB)--- = 
■ ■ ■ (PRX UB)(PR ZXZ ZUB)(PRX UB) ■■■ = ■■■ (PRX UB)(PR (-X) UB')(PRX UB) ■■■, which, by Theo- 
rem [5J is equivalent to A" 1 , since an overall phase factor of ±1 is irrelevant. Thus, the phase- flip channel is symmetric 
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FIG. 11: (color online) Onset of classicality is seen to be accentuated in a (Hadamard) quantum walk subjected to generalized 
amplitude damping with increasing temperatures. Figure (|9{a)) depicts the T — case (x — 1 in Eq. (|26[) ). (a) Finite 
temperature corresponding to x — 0.75 (b) T — oo, corresponding to x = 0-5- It may be noted that even at T — oo, for 
sufficiently small coupling the distribution remains non-classical. 




FIG. 12: (color online) Variation of standard deviation with amplitude damping noise level for various value of 9, 15°, 30°, 
45°, 60° and 75°. Note that the standard deviation for complementary angles converge to the same value. 



w.r.t. the operations Z and PRX. 

Regarding the bit-flip channel (|23|) : as in the above case, consider an unravelling X 2 = ■ ■ ■ (U B)(XU B)(U B) 
This is the same as: • • • (U B)(U* B")(U B) • • • , where, as may be seen by direct calculation, B" — XB. Now, ZX 2 = 
■■■(Z UB){X Z UB){Z UB) ■■■ = ■■■(Z UB)(X ZX XUB)(Z UB) ■■■ = ■ ■ ■ (Z UB)((-Z) WB"){Z UB) 
which, by Theorem [21 is equivalent to X 2 , since an overall phase factor of ±1 is irrelevant. 

Further, PRXl 2 = ••• (PRX UB)(X PRX UB){PRX UB) ■■■ = ... (PRX UB)(PRX X UB)(PRX UB) ■■■ = 
■ ■ ■ (PRX UB)(PRX WB")(PRX UB) • • • , which, by Thcoremd is equivalent to X 2 . ■ 



B. Decoherence via generalized amplitude damping channel 



Here we study the behaviour of quantum walk subjected to a generalized amplitude damping (with temperature 
T > 0), which would reduce at T = to the amplitude damping channel. As an example of a physical process 
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that realizes the generalized amplitude damping channel, we consider a two-level system interacting with a reservoir 
of harmonic oscillators, with the system-reservoir interaction being dissipative and of the weak Born-Markov type 
[2ll . I22I ] leading to a standard Lindblad equation, which in the interaction picture has the following form [2(| 



dt 1 



\t)=Y J {2R ] p s R]-R]R ] p s 



3=1 



(24) 



where i?i = (-y (N t h + l)/2) 1 / 2 R, R 2 = (7o^Vth/2) 1/2 i? t and N th = (exp(hu)/k B T) - l)" 1 , is the Planck distribution 
giving the number of thermal photons at the frequency ui, and 70 is the system-environment coupling constant. 
Here R = cr_ cosh(r) + e l ®a+ sinh(r), and the quantities r and are the environmental squeezing parameters and 
°"± = I ( a i i * cr 2)- For the generalized amplitude damping channel, we set r = <f> = 0. If T = 0, so that iVth = 0, 
then R 2 vanishes, and a single Lindblad operator suffices. 

The generalized amplitude damping channel is characterized by the following Kraus operators pol | : 



1 



y/1 -p(t) 



E 2 = ^l^ 



1 



E 3 = 




(25) 



where 



p(t) EE 1 - e -7o(2W th +l)t. 



1 



2^Vth + 1 



(26) 



When T = 0, x = 1, and for T -> 00, x = 1/2. 

The density operator at a future time can be obtained as [2C 



1(1 + A X ) A 2 



where 



^2 i(l-^i) 



A - e -7o(2Wth + l)t^ 



(27) 



(28) 



A 2 = e-^C^+D*^^))). (29) 

Figures [H [10] and [11] depict the onset of classicality with increasing coupling strength (related to p) and temperature 
(coming from x) ■ Figure [5] which shows the effect of an amplitude damping channel on a Hadamard walker at zero 
temperature, illustrates the breakdown of RX symmetry even though the initial state is |0) + i|l). This is because, 
in contrast to the phase-flip and bit-flip channels, the generalized amplitude damping is not symmetric towards the 
states |0) and |1). However, the extended symmetry PRX is preserved both for Hadamard as well as biased walks, 
as seen from Figures [9] and [TO] respectively. 

From Figures ©a)), l)ll( a.b)). the onset of classicality with increasing temperature is clearly seen. Figure [T2l 
presents the standard deviation for quantum walks on a line with various biases, subjected to amplitude damping 
noise. The standard deviation for complementary angles (8 <-> ir/2 — 9) is seen to converge to the same value in the fully 
classical limit. This may be understood as follows. First, wc note that since PRX is a symmetry of the quantum walk, 
and the effect of P does not show up in the standard deviation plots, RX by itself is an apparent symmetry. Further, in 
the classical limit the measurement outcome being a unique asymptotic state for the (generalized) amplitude damping 
channel, effectively X ~ 1, which makes R a symmetry operation. 

The following theorem generalizes Theorem [5] to an open system subjected to a generalized amplitude damping 
channel. 

Theorem 5 The operations Z and PRX are symmetries for the generalized amplitude damping channel. 

Proof. By virtue of Theorem[3] it suffices to show that any given unravelling is invariant under Z and PRX. Consider 
an unravelling 



X 3 = ■■■(E UB)(E 1 UB)(E 2 UB)(E 3 UB)--- 
= ■■■(UB {0) )(U^B {1) )(UB {2) )(U^B {3) )--- , 



(30a) 
(30b) 
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where the non- unitary matrices are given by B(j\ — EjB. Now, 

ZX 3 = ■■■(E a ZUB)(E 1 ZUB)(E 2 ZUB)(E 3 ZUB)--- 

= ■■■ {ZE UB)((-Z)E 1 UB)(ZE 2 UB)((-Z)E 3 UB) ■ ■ ■ 

= ■■■ (ZUB (0) )((-Z)U^B (1) )(ZUB (2) )((-Z)U^B {3) ) ■ ■ ■ 

= ■■■ (UB^i-U^B^iUB^X-U^B^) ■ ■ ■ (31) 

Ignoring the overall ±1 factor in Eq. ([3"T]) . and comparing it with Eq. (|30a[) . and noting that that the derivation of the 
proof of Theorem [2] did not require the matrices Bj to be unitary, we find along similar lines that ZX 3 is equivalent 
to X 3 . 

The following may be directly verified 

PRX1 3 = • • • (E PRXUB)(E 1 PRXUB)(E 2 PRXUB)(E 3 PRXUB) ■ ■ ■ (32a) 
= ■■■{E UB (2) ){E l UB {2) ){E 2 UB {2) ){E 3 UB (2) )--- (32b) 
= • •■ (UB i2 })(U^B^)(UB (2 })(U^B^) • • • , (32c) 

which, by Theorem^ is equivalent to X 3 . For proof of Eq. (|32b|) , see the proof of Theorem [TJ Eq. (|32cp is obtained 
analogously to Eq. (|30b[) . except that the matrix B^> is used instead of B. ■ 

This may be expressed by the statement 

NW ~ AfZW, (33a) 
AfW ~ AfPKXW, (33b) 

which generalizes Eq. (|10p . These results show that the symmetries persist for dephasing (phase flip), bit flip and 
(generalized) amplitude damping channels. 

V. QUANTUM WALK ON A CYCLE 

In this work, though we are primarily concerned with symmetries for a quantum walk on a line, and the influence of 
noise on them, we briefly consider in this Section an extension of the above ideas to quantum walks on a cycle. Further 
extensions would be quantum walks on a more general graph [23l l2~ij or in higher dimensions d > 2 [25[ . In the former, 
the ID walk is generalized to TV-cycles and to hypercubes, including the effect of phase noise in the coin space, and 
decoherence in position space. In the latter, the Hadamard transformation is generalized to a non-entangling tensor 
product of Hadamards, or to an entangling discrete Fourier transform or the Grover operator. They bring in many 
novel features absent in the quantum walk on a line. Here we will restrict ourselves to pointing out that quantum 
walk on a cycle differs considerably from walk on a line, both with respect to symmetry operations as well as noise. 

In contrast to the case of quantum walk on a line, none of the four discrete symmetries of Theorem Q] hold in general 
for unitary quantum walk on a cycle or closed path. Thus, if B in Eq. © is replaced by any of B^\ B^ 2 \ B^ 3 \ or 
B^\ given by Eq. ([8]), the spatial probability distribution is not guaranteed to be the same. 

Theorem 6 The operation G : B — > B* is in general not a symmetry of the quantum walk on a cycle. 
Proof. For the cyclic case, in place of Eq. ©, we now have 

= (UB) n \a,/3)= Yl K,j n -i---bh,hb jua \j n ,0 + 2J-n(modR)), (34a) 

31, ht— Jn 

|* a ) = (UB^) n \a,(3)= £ b 3n ^_ x ■ ■ ■ b h , h b jl , a (e^y^+-+^+°'\j n , (3 + 2J - n (mod R)), (34b) 

jlda,— ,jn 

where R is the number of sites in the cycle. For an arbitrary state \a,b) in the computational-and-position basis, we 
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FIG. 13: (color online) An instance of breakdown of phase flip symmetry in a unitary quantum walk on a cycle, where the 
two extreme points (located at positions ±50) on the plot are spatially adjacent. The number of sites is 101 and t = 5000 (in 
units of discrete time-steps), with bias angle 9 = 30°. (a) The solid curve represents the positional probability distribution 
without any symmetry operation applied, while the dashed curve represents that with a phase flip operation applied at each 
walk step, (b) The same as the above, but with time-averaging applied over every 50 steps, in order to more clearly bring out 
the breakdown in symmetry. 



have 



(a,6|*i) = ^ ,K '-J- : ■■■t'j.'J ( 35a ) 

jl,j2,— ,jn-l£J 

E Ko n -,---K,hKA^% (35b) 



where J is the set of binary (n — l)-tuples j±, ]2, ■ ■ ■ ,j n -i such that J = j\ + j% + • • • + j n -i + a satisfies b = 
(3 + 2 J - n (mod R). We find that e = a - a + J and J = (b + n - (3)/2 + mR, where m = 0, 1,2,- • • , [n/R\ . 
Thus, the terms in the superposition (|35ap are not in general identical with those in (|35b[) . apart from a common 
factor, unless = 0, 2ir, Air, • • • . A similar argument can be used to show that the terms in (a, 6|\E'i) are not in general 
the same as those in (a,b\^2)- Given the independence of <j> from the coefficients bj 2 j 1 bj 1>a ^ it is not necessary that 
|(a, 6|\&i)| 2 + |(o, bl^i)! 2 = |(a,6|*2)| 2 + |(a,&|*2)| 2 - The equality holds in general (for arbitrary unitary matrix B 
and time parameter n) if and only if 4> = 0, 2tt, 4n, ■ ■ ■ . Repeating the argument for B^ 2 \ and B^ 4 \ we find that 
all the four discrete symmetries of Theorem [T] break down in general. ■ 

We note that for phase flip symmetry, where e 1 ^ — — 1, the superposition terms in Eq. (|35b[) may differ from the 
corresponding terms in Eq. (|35aP only with respect to sign. Given that all the terms like ij\ 5 ' cfcCi are built 
from a small set of trignometric functions of the three parameters 8, ( and £, certain values of n may render the right 
hand sides of Eqs. ()35a|l and (I35bj) equal. However in general, this equality will not hold for arbitrary n. 

An instance of breakdown of phase flip symmetry in the unitary quantum walk on a cycle is demonstrated in the 
example of Figure 1131 The profile of the position probability distribution varies depending on the number of sites 
and the evolution time. Remarkably, this symmetry is restored above a threshold value of noise. The pattern in 
Figure [14] corresponds to phase noise with p = 0.02 applied to a quantum walk, either with or without a phase flip 
symmetry operation. We note that the introduction of noise tends to classicalize the random walk, hence causing 
it to asymptotically reach a uniform distribution (23l . [24| • The above-mentioned symmetry restoration happens well 
before the uniformity sets in. The initial lack of symmetry gradually transitions to full symmetry as the noise level is 
increased. Thus, the role of symmetry operations and noise is quite different in the case of quantum walk on cycles 
as compared with that on a line. 
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FIG. 14: (color online) Restoration of phase flip symmetry in a noisy quantum walk on a cycle, where the two extreme points 
on the plot are spatially adjacent. The number of sites is 101 and t = 5000 (in units of discrete time-steps), with bias angle 
6 = 30° . The figure depicts the position probability distribution with or without phase flip symmetry operation applied at each 
step, with phase damping noise level p = 0.02 [Eq. (|18[) ]. After sufficiently long time, the quantum walk reaches the uniform 
distribution, typical of classical random walk. 



A more detailed treatment of symmetries and noise in a quantum walk on a cycle and general graphs of other 
topologies will be presented elsewhere [26| . 

VI. EXPERIMENTAL REALIZATION IN PHYSICAL SYSTEMS 

Experimental realization of quantum walk using any of the proposed schemes is not free from noise due to envi- 
ronmental conditions and instrumental interference. In particular, noise can be a major issue in the scaling up of the 
number of steps in already realized quantum walk systems. Understanding the symmetries of the noisy and noiseless 
quantum walk could greatly help in the improvement of known techniques and in further exploration of other possible 
systems where quantum walk can be realized on a large scale. In this section we discuss the realization in a nuclear 
magnetic resonance (NMR) quantum-information processor and ultra-cold atomic systems. 



A. NMR quantum-information processor 

Continuous time [Tl| and discrete time [Toj quantum walk have been successfully implemented in a nuclear mag- 
netic resonance (NMR) quantum-information processor. Considering the benefits of the effect of decoherence on the 
quantum walk [27], Hfl , Ref. [l(| has also experimentally added decoherence on the discrete time quantum walk by 
implementing dephasing in NMR. By adding decoherence after each step they have shown the transition of quantum 
walk to the classical random walk. 

In NMR spectroscopy of the given system (molecule) the extent of its isolation from the environment is determined 
in terms of its phase coherence time and its energy relaxation time T\. If the pulse sequence is applied to the NMR 
quantum- information processor within the time T < T^,Tx, the system is free from the environmental effects. The 
pulse sequence exceeding the time Ti can be considered to affected by the dephasing channel and the pulse sequence 
exceeding the time T\ can be considered to be affected by the amplitude damping channel. In experiments of time 
scale greater than the time Ti or T\, a refocusing pulse sequence is applied to compensate for the environmental 
effects. In Ref. [Tcj the pulse sequence for the quantum walk was implemented within the time T\ and T<i. 

The environmental effect (noise) on quantum walk symmetries presented in this article can be verified in the NMR 
system by scaling up the number of steps of quantum walk realized. By applying a controlled amount of the refocusing 
pulse sequence, the effect of different levels of noise can be experimentally verified. 
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B. Ultra-cold atoms 

There have been various schemes suggested to implement quantum walks using neutral atoms in an optical lattice 
EH • In Ref . , the controlled coherent transport and splitting of atomic wave packets in spin dependent optical 
lattice has been experimentally demonstrated using rubidium atoms. A Bose-Einstein condensate of up to 3 x 10 5 
atoms is initially created in a harmonic magnetic trap. A three dimensional optical lattice is superimposed on the 
Bose-Einstein condensate and the intensity is raised in order to drive the system into a Mott insulating phase [30l ] . 

Two of the three orthogonal standing wave light fields is operated at one wavelenght, X VyZ = 840 nm and the third 
along the horizontal direction is tuned to the wavelength X x — 785 nm between the fine structure splitting of the 
rubidium Dl and D2 transitions. Along this axis a quarter wave plate and an electro-optical modulator (EOM) is 
placed to allow the dynamic rotation of the polarization vector of the retro-reflected laser beam through an angle 9 by 
applying an appropriate voltage to the EOM. After reaching the Mott insulating phase, the harmonic magnetic field 
is completely turned off but a homogeneous magnetic field along the x direction is maintained to preserve the spin 
polarization of the atoms. The light field in the y and z direction is adiabatically turned off to reduce the interaction 
energy, which strongly depends on the confinement of the atoms at a single lattice site. 

A standing wave configuration in the x direction is used to transport the atoms. By changing the linear polarization 
vector enclosing angle 9, the separation between the two potentials is controlled. By rotating the polarization angle 
9 by 7r, with the atom in a superposition of internal states, the spatial wave packets of the atom in the |0) and the 
|1) state are transported in opposite directions. The final state after such a movement is then given by l/v^dOjX — 
1) + iexp(i/3 p )|l, x + 1)). The phase (3 P between the separated wave-packets depends on the accumulated kinetic and 
potential energy phases in the transport process and in general will be nonzero. The coherence between the two 
wave-packets is revealed by absorption imaging of the momentum distribution. A 7r/2 microwave pulse is applied 
before absorption imaging to erase the which-way information encoded in the hyperfine states. 

However, to increase the separation between the two wave-packets further, one could increase the polarization angle 
9 to integer multiples of 180°. To overcome the limitation of the maximum voltage that can be applied to the EOM, 
a 7r pulse after the polarization is applied, thereby swapping the roles of the two hyperfine states. The single particle 
phase (3 P remains constant throughout the atomic cloud and is reproducible. After the absorption imaging a Gaussian 
envelope of the interference pattern is obtained. 

One can build up on the above technique to implement a quantum walk which introduces phase along with each 
splitting (step). The above setup can be modified by dividing the separations into small steps and introducing a 7r/2 
pulse after each separation without intermediate imaging. A phase (3 P is introduced in each step. The absorption 
imaging of the distribution of the atomic cloud after n steps would give the interference pattern of the quantum walk. 

This effect of the addition of phase during the quantum walk process can be easily understood from the phase 
damping channel and arbitrary phase rotation presented in this article. The addition of 7r pulse to overcome the limit 
of EOM is a bit flip operation in the quantum walk. 

There have been other proposals for physical realization of quantum walk using Bose-Einstein condensate (BEC) 
[l5T | where the unitary shift operator induces a bit flip. A stimulated Raman kick is used as a unitary shift operator 
to translate the Bose-Einstein condensate in the Schrodinger cat state to a superposition in position space. Two 
selected levels of the atoms in the Bose-Einstein condensate are coupled to the two modes of counter-propagating 
laser beams. The stimulated Raman kick, in imparting a translation in position space, also flips the internal state of 
the Bose-Einstein condensate. An rf pulse (tt pulse) is suggested as a compensatory mechanism to flip the internal 
state of the condensate back to its initial value after every unitary shift operator. From the PRX symmetry pointed 
in this article, it follows that there is no need for the compensatory operation for an unbiased quantum walk started 
in the state |0) + i\l). The availability of walk symmetries could also be useful for exploring other possible physical 
implementations which induce such symmetry operations along with the translation. 

In the most widely studied version of quantum walk, a quantum coin toss (Hadamard operation) is used after every 
displacement operation. Continuous external operations on a particle confined in a trap reduces the confinement time 
of the particle. Reducing the number of external operations will benefit the scaling up of the number of steps of 
the quantum walk. One can workout a system where the quantum coin toss operation is eliminated by transferring 
the burden of evolving the particle in superposition of the internal state to the displacement operator itself [3lL I32I ] . 
Such a transfer may introduce additional operations that correspond to the walk symmetry, and may thus be ignored. 
Further, systems of this kind are expected to be affected by amplitude damping as one of the states might be more 
stable than the other one in the trap. The present work could help to optimize such a noisy quantum walk. 

In a scheme suggested using quantum accelerator mode fl6j . different internal states of an atom receive different 
momentum transfer with each alternative kick, giving different walking speeds in the two directions and this can 
be seen as a biased walk. Our study of symmetry in a noisy, biased quantum walk could help improve the above 
technique and make it easier for its experimental realization. 
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VII. CONCLUSIONS 

Our work considers variants of quantum walks on a line which are equivalent in the sense that the final posi- 
tional probability distribution remains the same in each variant. In particular, we consider variants obtained by the 
experimentally relevant operations of Z or X applied at each quantum walk step, with the symmetry operations 
given by Z and PRX. This could be experimentally advantageous since practical constraints may mean that one 
of the variants is preferred over the rest. A specific example is the simplification of the implementation of quantum 
walk using a Bose-Einstein condensate with a stimulated Raman kick providing the conditional translation operation. 
What is especially interesting is that these symmetries are preserved even in the presence of noise, in particular, those 
characterized by the phase flip, bit flip and generalized amplitude damping channels. This is important because it 
means that the equivalence of these variants is not affected by the presence of noise, which would be inevitable in 
actual experiments. The symmetry of the phase operation under phase noise is intuitive, considering that this noise 
has a Kraus representation consisting of operations that are symmetries of the noiseless quantum walk. However, for 
the PRX symmetry under phase noise, and for any symmetry under other noisy channels (especially in the case of 
generalized amplitude damping channel), the connection was not obvious before the analysis was completed. 

Our results are supported by several numerical examples obtained by evolving the density operator in the Kraus 
representation. However, analytical proofs of the effect of noise on symmetries are obtained using the quantum 
trajectories approach, which we find convenient for this situation. We also present the quantum walk on the cycle, 
which can be generalized to any closed graph. An interesting fact that comes out is that the symmetry breaks down 
in general but is restored above a certain noise level. Some representative plots demonstrating the effect of the phase 
damping channel on phase flip symmetry are presented. 

Finally we have discussed the experimental realization of quantum walk in a physical system, such as NMR and 
ultra cold atoms, as examples where these studies can be beneficial for improving the efficiency of the implementation 
on large scales. 
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